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Abstract 

The purpose of this paper is to provide a brief review of some recent developments in quantum feedback 
networks and control. A quantum feedback network (QFN) is an interconnected system consisting of open 
quantum systems linked by free fields and/or direct physical couplings. Basic network constructs, including 
series connections as well as feedback loops, arc discussed. The quantum feedback network theory provides 
a natural framework for analysis and design. Basic properties such as dissipation, stability, passivity and 
gain of open quantum systems are discussed. Control system design is also discussed, primarily in the 
context of open linear quantum stochastic systems. The issue of physical realizability is discussed, and 
explicit criteria for stability, positive real lemma, and bounded real lemma are presented. Finally for 
linear quantum systems, coherent H°° and LQG control are described. 

Key Words Open quantum systems; quantum feedback networks; physical realizability; H x control; 
LQG control. 

1 Introduction 

Quantum technology is an interdisciplinary field that studies how to engineer devices by exploiting their 
quantum features. Regarded as the second quantum revolution, quantum technology has many potential 
far-reaching applications [10]. For example, Shor [IB] presented a quantum algorithm which can offer 
exponential speedup over classical algorithms for factoring large integers into prime numbers. Bennett et 
al. [3] proposed a quantum teleportation protocol where an unknown quantum state can be disembodiedly 
transported to a desired receiver. Atomic lasers hold promising applications in nanotechnology such as 
atom lithography, atom optics and precision measurement [33]. Quantum technology (including quantum 
information technology) has more powerful capability than traditional technology and is one of the main 
focuses of scientists. Nevertheless, many challenging problems require to be systematically presented and 
successfully addressed in order to foster wider real- world applications of quantum technology in our life [10] . 

Recent years have seen a rapid growth of quantum feedback control theory [301 HJ [9l G2 l24| 160] , If 
measurement is involved in the feedback loop, the feedback mechanism is conventionally called measurement- 
based feedback, e.g., [2j [581 El El E E21 [351 ED]. Measurement-based feedback control of quantum systems is 
important in a number of areas of quantum technology, including quantum optical systems, nano-mechanical 
systems, and circuit QED systems. In measurement-based feedback control, the plant is a quantum system, 
while the controller is a classical (namely non-quantum) system. The classical controller processes the 
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outcomes of measurement of an observable of the quantum system (e.g. the number of photons of an optical 
field) to determine the classical control actions (e.g. magnetic field) that are applied to control the behavior 
of the quantum system. Classical controllers are typically implemented using standard analog or digital 
electronics. However, for quantum systems that have bandwidth much higher than that of conventional 
electronics, an important practical issue for the implementation of measurement-based feedback control 
systems is the relatively slow speed of standard classical electronics, since the feedback system will not work 
properly unless the controller is fast enough. 

Alternatively, quantum components may be connected to each other without any measurement devices 
in the interconnections. For example, two optical cavities can be connected via electromagnetic fields (light 
beams). Such feedback mechanism is referred to as coherent feedback as originally proposed in [65 1 159 1 [36 ], 
[2"71 loTl E2] ■ The interconnection of a quantum plant and a quantum controller produces a fully quantum 
system; quantum information flows in this coherent feedback network, thus coherence is preserved in the 
whole quantum network. Moreover, a coherent feedback controller may have the similar time scale as the 
plant, and likely would be much faster than classical signal processing. Finally, it is becoming feasible 
to implement quantum networks in semiconductor materials, for example, photonic crystals are periodic 
optical nanostructures that are designed to affect the motion of photons in a similar way that periodicity 
of a semiconductor crystal affects the motion of electrons, and it may be desirable to implement control 
networks on the same chip (rather than interfacing to a separate system) , |33[ l43l [49l l44"l 150] . For several 
reasons, then, it is desirable to implement controllers using the same or similar (e.g. in time scales) hardware 
as plants. Therefore, it might be advantageous to design coherent feedback networks. 

Recently there is a growing interest in the study of coherent quantum feedback networks and control. 
For example, quantum feedback network structure has been studied [591 ED CES1 El [231 CE1 ES] . The prob- 
lem of H°° control has been discussed [221 [3"1~1 [28l 166] . The problem of coherent LQG control has been 
investigated [371 EH]- The issue of physical realizability has been analyzed [221 14*6] [29l I66| . The problem of 
network synthesis of quantum systems via optical devices have been studied |38[ 139 1 1^1 142j . There are also 
many papers investigating the applications of coherent feedback control, such as intra-cavity squeezing [59], 
optical field squeezing [62j "TTJ [2Tj , H°° control [31], entanglement enhancement [19], polarization squeezing 
[471 [45] , error correction in quantum memories [251 [26] , and optical switches [32] . 

The paper is organized as follows. Section [2] discusses closed quantum systems, in particular, closed 
quantum harmonic oscillators. Section[3]introduces Boson field and a basic model structure of open quantum 
systems. Section 0] discusses mechanisms by which open quantum systems interconnect. Section [6] focuses 
on linear quantum systems. Section [5] presents results for quantum dissipative systems. Section [7J discusses 
fundamental characteristics of linear quantum systems. Section [8] presents H°° and LQG controller synthesis 
of linear quantum systems. Section [9] touches on how linear quantum systems can be realized by means of 
optical devices. Section [10] concludes the paper. 

Notation, i is the imaginary unit. 5jk is Kronecker delta, and 8(t) is Dirac delta. Given a column 
vector of operators or complex numbers x = [ x\ ■ ■ ■ x m ] T where m is a positive integer, define x& = 
* ] ) where the asterisk * indicates Hilbert space adjoint or complex conjugation. Denote x* = 
(x#) T = [ x* ■■■ x* m ]. Furthermore, define the doubled-up column vector to be x = [ x T (x^) T ] T ■ 
The matrix case can be defined analogously. Given two matrices U, V £ C rxfc , a doubled-up matrix A (U, V) 
is defined as A(U,V) := [ U V; ]. Let I n be an identity matrix. Define J n = diag(/ n ,— I n ) 

and @ n = [0 I n ; — I n 0]. (The subscript "n" is always omitted.) Then for a matrix X £ C 2nx2m , define 

:= J m X^ J n . Finally we also use / to denote identity operators. 
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2 Closed Systems 



In this paper, closed systems means systems that have no interactions with other systems and/or environ- 
ment. In this section starting from the fundamental Schrodinger's equation for closed quantum systems, we 
introduce closed quantum harmonic oscillators which in later sections will be allowed to interact with other 
systems or electromagnetic fields to produce open quantum systems. 

Given a closed quantum system with Hamiltonian H, we have the following Schrodinger's equation 

^U(t) = -iHU(t), 17(0) = I. (1) 

Clearly, U(t) is a unitary operator. The system variables X(t) evolve according to X(t) = U*(t)XU(t) with 
initial point X(0) = X, which satisfy, the Heisenberg picture, 

±X{t) = -i[X(t),H(t)]. (2) 

(Note that for closed systems H(t) = H for all t due to preservation of energy). 

Alternatively, the system density operator p(t) = U(t)pU*(t) with p(0) = p satisfies, the Schrodinger 
picture, 

±p(t) = -i[H,p(t)}. (3) 
2.1 Closed Quantum Harmonic Oscillators 




Figure 1: Closed optical cavity, black rectangles denote fully reflecting mirrors at cavity resonant frequency 

An example of closed quantum harmonic oscillators is an optical cavity with H = uja*a (upon scaling), 
Fig. [TJ where uj is the resonant frequency, and the annihilation operator a is the cavity mode (an operator on 
a Hilbert space). The adjoint operator a* of a is called the creation operator, a and a* satisfy the canonical 
commutation relation [a(t),a*(t)] = 1 for all t > 0. Finally by Eq. (J2J), 

-^■a(i) = -iua(t), a(0) = a. (4) 

Therefore a(t) = e~ lujt a — an oscillator. 

In general, let G be a closed quantum system of interconnection of n quantum harmonic oscillators. The 
behavior of G is determined by the Hamiltonian 

a, (5) 

where 0_ and 0+ are respectively C nxn matrices satisfying f2_ = Q\_ and £7+ = fiT. By Eq. ([2]), 

aj(t) = -i[aj(t),H (t)}, Oj-(0) = aj, (j = l,...,re). 
lr The reduced Planck constant h is omitted throughout the paper. 



i7n 



n* n* 



3 



In a compact form we have the following linear differential equations 

a(t) = A Q a(t) 

with initial condition a(0) = a, where 

A = -A(ifi_,ifi+). 



(6) 
(7) 



3 Quantum Fields and Open Quantum Systems 
3.1 Boson Fields 

The m-channel Boson field b(t) = [bi(t), . . . , b m (t)] T are operators on a Fock space T [U], whose components 
satisfy the singular commutation relations 



[b j (t),b* k (t')]=6 jk 6(t-t'), [b 1 (t),b k (t')]=0, [b*(t),bl(t')}=0, (j,k = l,...,m). 



(8) 



The operators bj(t) may be regarded as quantum stochastic processes, see, eg., [131 Chapter 5]; when the field 
is in the vacuum state, namely absolutely zero temperature and completely dark, they are called standard 
quantum white noise (that is, M = N = in \13\ Eq. (10.2.38)]). The integrated processes Bj(t) = 
Jo ^i( r )^ r are quantum Wiener processes with Ito increments dBj{t) = Bj(t + dt) — Bj(t), (j = 1, . . . , m). 
There might exist scattering between channels, which is modeled by the gauge process 



A(t)= / b*(T)b T (T)dl 



An(t) 



A m i(i) 



A im (t) 



Amm (t) 



(9) 



with operator entries Aj k on the Fock space T . Finally in this paper it is assumed that these quantum 
stochastic processes are canonical, that is, they have the following non-zero Ito products 



dBj(t)dB* k (t) = Sjkdt, dA jk dBi(t) = 5 k idB*(t), 

dBj(t)dA kl (t) = 5 jk dBi(t), dA jk {t)dAi m {t) = 6 k idA jm (t), (j,k,l = 1, ...,m) 
3.2 Open Quantum Systems in the (S, L, H) Parametrization 



(10) 



Figure 2: Open optical cavity, the white rectangle denotes a partially transmitting mirror at cavity resonant 
frequency 

When a quantum system G is driven by a Boson field J 7 , we have an open quantum system. For 
example, if we allow the closed optical cavity in Fig. [1] to interact with a Boson field, we end up with an 
open optical cavity, Fig. [2l While the mutual influence between the system and field may be described 
rigorously from first principles in terms of an interaction Hamiltonian, it is much more convenient to use 
an idealized quantum noise model which is valid under suitable rotating wave and Markovian assumptions, 
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as in many situations in quantum optics, eg., cascaded open systems, see 021 El EI] for detail. Let Aq 
and Af be physical variable spaces of the system G and the field J- respectively, then the physical variable 
space for the composite system is the tensor product space Aq <8> Ap- 

Open quantum systems G studied in this paper can be parameterized by a triple (S, L, H) |2Q|. I16J . Here, 
S is a scattering matrix with entries in the system space Ac, L € Ac is an coupling operator that provides 
interface between systems and fields, H £ Ag is the internal Hamiltonian of quantum system G. 

With these parameters, and assuming that the input field is canonical, that is, Eq. (|10j) holds, we have 
the following Schrodinger's equation for open quantum systems (in Ito form) 

dU(t) = |tr[(5 - I m )d\ T ] + dB\t)L - L ] SdB{t) - {iH + ^L)dt\ U(t), U(0) = I. (11) 

(Note that for a closed system, Eq. (jlip becomes the familiar Schrodinger's equation ([1]).) This, together 
with the evolution X(t) = U(t)*XU(t), yields the following quantum stochastic differential equations (QS- 
DEs), in Ito form, 

dX{t) = {-%[X{t),H{t)\+C L{t) {X{t)))dt+ (12) 
+tv[{S\t)X(t)S{t) - X(t))dA T (t)}, X(0) = X, 

where the Lindblad operator Cl is 

C L (X):=h^[X,L] + ^,X]L. (13) 

For later use, we define a generator operator 

g G {X):=-i[X,H]+C L (X). (14) 

The output field B out (t) = U* (t) B (t)U (t) satisfies 

dB out (t) = L(t)dt + S(t)dB(t). (15) 

The gauge process of the output field A otti (t) := /* 6* t (s)6^ ut (s)ds = U*(t)A(t)U(t) satisfies 

dA out (t) = S#(t)dA(t)S T (t) + S#(t)dB*(t)L T (t) + L*(t)dB T (t)S T (t) + L*(t)L T \t)dt. (16) 

Finally, in the Schrodinger picture, the reduced system density operator p satisfies the master equation, cf. 
[13 Sec. 11.2.5], 

j t p{t) = -i\H,p{t)]+L' L {pm (17) 
where the operator C' L is defined to be 

C' L (p) := L T pL# - h)Lp - l -pL^L. (18) 
Remark 1 Clearly, open quantum systems presented in this section are quantum Markov processes. 



5 



3.3 Examples 
3.3.1 Optical Cavity 

The one degree of freedom closed quantum harmonic oscillator in Fig. Q]can be described by (— , - ,wa*a), 
where the symbol "-" means that there is neither scattering nor coupling. The open optical cavity in 
Fig. [5] may be described by (1, ^/Jta, uo*o), where k is coupling coefficient and u) is the resonant frequency. 
According to Eqs. (fl~2>(fT5|h 



da(t) = 
dB 0Ut (t) = y/na(t)dt + dB(t) 



iuj + —)a(t)dt — \/7idB(t), a(0) = a, 



3.3.2 Two-level Systems 

Given a two-level system parameterized by 



S- = 1, L = y/Kcr-, H = -a z , 



with Pauli matrices 



a. 
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by Eq. (fTT|) the unitary operator U(t) evolves according to 

dU(t) = {^dB*(t)a^ - ^<r + dB{t) - j(a z + l)dt - i^a z dt} U(t), U(0) = I. 
By Eqs. O and fTH]l 

da x (t) = -(^a x {t) + ua v (t))dt + V^dB*(t)a z {t) + s/na z {t)dB{t), 



da y {t) = —(—a y (t) — uja x (t))dt — iy/KdB*(t)a z (t) + i\/ r Ka z (t)dB(t), 



da z (t) 



-k(I + a z (t))dt 



'o x (t)-io y {t))dB*{t) 



(a x (t) + iay(t))dB(t). 



On the other hand, the output field is 



dB out {t) = y/K<T-(t)dt + dB(t), 



(19) 
(20) 



(21) 

(22) 

(23) 
(24) 

(25) 

(26) 
(27) 



dA out (t) = dA(t) + y/KdB*{t)a-(t) + y/~K<j + {t)dB{t) + ~(cr z + l)dt. 
Remark 2 It can be seen from Eq. ([23]) - ([25]) that two-level systems are nonlinear quantum systems. 

4 Interconnection 



In this section we discuss how two quantum systems can be connected to each other. More specifically 
we discuss concatenation product, series product, direct coupling, and linear fractional transform. Several 
examples from the literature are used to illustrate these interconnections. Propagation delays are ignored 
in interconnections. Discussions of influence of propagation delays on system performance can be found in, 
eg., [21]. 
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Figure 3: Concatenation product G\ EE G% 



Figure 4: Series product G2 < G\ 



4.1 Concatenation Product 

Given two open quantum systems G\ = (S±, L\, Hi) and G2 = (S2, L2, H2), their concatenation product, 
Fig. [3j is defined to be 



G\ EH G2 '■- 





" 5i 







" £1 " 


( 





s 2 


J 


£2 



Hi + H 2 



(28) 



4.2 Series Product 

Given two open quantum systems G\ = (S±, L%, Hi) and G2 = (5*2, L2, H2) with the same number of input, 
their series product, Fig. U is defined to be 



G 2 < d := \ S2Su L 2 + StLuH! + H 2 + ^(L^Lx - L\slL 2 )j . 



(29) 



Theorem 1 (Principle of Series Connections, ]16\ Theorem 5.5]): The parameters of the composite system 
(?2 <— G\, obtained from G1EHG2 when the output of G\ is used as input ofG 2 , is given by the series product 
G 2 <Gl 

4.3 Direct Coupling 

In quantum mechanics, two independent systems G\ and G2 may interact by exchanging energy. This energy 
exchange may be described by an interaction Hamiltonian Hi n t of the form Hi n t = x\X2 + X^Xi, where 
X\ G Agy an d X 2 G -A.G21 see i e g-> EH]) [2Z], [66]. In this case, we say the two systems G\ and G2 are 
directly coupled, and the composite system is denoted G\ \x\G2, Fig. [5j 



4.4 Linear Fractional Transform 

Let G in Fig. [6] be of the form 



S\\ S\2 

S21 S22 



Li 
L 2 



,H \. Assume that (I — S22) 1 exists. Then the 



conventional linear fractional transform yields a feedback network 

F(G) = (S n +Si2(I-S22y 1 S2i,L 1 +S 1 2(I-S22y 1 L 2 ,H+Im{L{S 1 2(I-S22y 1 L2}+lm{LlS22(I-S22y 1 L2}) 
from input b to output b ou t- 
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Figure 5: Directly coupled system G\ txi G2 

b 



Figure 6: Linear fractional transformation F(G) 



4.5 Examples 

In this section examples in the literature are used to demonstrate the usefulness of the parametrization 
(S,L,H) and interconnections. More examples can be found in |59 | 1 16 ^ HH [ I32j. For the convenience of the 
readers to refer to the original papers we use symbols in those original papers. 

Example 1 ([6]) 

In [6] quantum trajectory theory is formulated for interaction of open quantum systems via series product, 
Fig. [H Given two open systems G± = (1, La, Ha), G2 = (1, Lb, Hb) with La = \J2kao>a, Lb = V^bo-b- 
Here ka and kb are coupling constants and a a and as are annihilation operators for systems G\ and G2 
respectively. The series product yields 



G 2 < G\ = (1, \/2kacla + \JlK B a B , H A + H B + i\J k>akb ( a A a B ~ a* B aA))- 



(30) 



Identifying Ha + Hb + iy/KAKB~(a*A a B — a*B a A)) with Hs in [6j Eq.(7)] and ^/2ka~cla + V^bo-b with C in 
[61 Eq.(9)] respectively, Eq. (fT7|) re-produces the master equation [6l Eq.(8)]. 

Example 2 ([J2]) 



In [12] quantum Langevin equations and a quantum master equation were derived for a cascade of two 
two-level systems, Fig. [71 The two two-level systems Gi and G2 in Fig. [7J are respectively 

Gi = (1, y^r, 0) ffl ((1, y^Tar, 0) < (1, E, 0)), 

G 2 = (1 , V%^2 > 0) ffl (1 , Vt^ct 2 " , 0) . 

Here, 71 and 771 are coupling constants for system Gi, is an incident coherent electric drive (not an 
operator) of G\. 71 and rji are coupling constants for system G2. cr~ is defined in Eq. ([2T|) : 



0"n 
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Gi 
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Figure 7: A re-draw of [12j FIG.l]. For clarity, coupling operators • v /7i<t 1 , y^TTfcr-L , y/jia 2 , yfnicr-^ are shown 
explicitly. E is a coherent electric drive of G±. 



Let G be a series product G = (G 2 ffl (1,0,0)) < ((1,0,0) ffl Gi). With these, according to Eq. (UJ, [321 
Eq.(14)] can be re-produced. (Notice the fact that in the interaction picture H sys = is used.) 

Example 3 ( [471 135] ) 

A scheme is proposed in [47] to produce continuous-wave fields or pulses of polarization-squeezed light by 
passing classical, linearly polarized laser light through an atomic sample twice; that is, the field output of 
the first pass is fed back to the atomic sample again so as to generate polarization-squeezed light. This 
scheme is confirmed and extended in [35]. The atomic sample can be modeled as an open quantum system 
G = Gi ffl G 2 with 

1 i 
Gi = (1, -j=ap,0), G 2 = (l,--j=ax,0). 

Here a is a real constant, and x,p are position and momentum operators respectively, cf. |45] Eqs. (Al)- 
(A2)]. Double-pass of an electromagnetic field through the atomic sample introduces a series product, that 
is the overall quantum system is 

a c? 
G 2 < Gi = (I, -(p - ix), — (xp + px)), 

whose Schrodinger equation is (33 Eq.(l)]. 
Example 4 ((66]) 

Given two closed quantum Harmonic oscillators G\ and G 2 as studied in Section 12.11 we take the interaction 
Hamiltonian Hi n t in Fig. [5] to be 

^ = ^(a (1) W)+a( 2 )t-aW), (31) 

where S = A(iK-,iK+) for matrices K_,K + £ ^n 2 xn 1 ^ -p^g Hamiltonian for the directly coupled system 
Gi cxi G 2 is 

H = Hq\ + Hi n t + HQ y2 , (32) 

where ff ,fc = A(^\ nj } ) is the self-Hamiltonian for Gj~, and Hi nt is given by ([3T|) . It is easy 

to show that system operators a^\t) = U*(t)a^U(t) (j = 1,2) satisfy the following linear differential 
equations, in Stratonovich form, 

&M(t) = A),ia (1) W + M (2) (i), a«(0) = a«, 
ti 2 \t) = A , 2 a^(t) + B 21 a^(t), a^(0) = a^ 2 \ 
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Figure 8: Plant-exosystem network PAW 



where 

A j = -A(in^\in { i } ), B l2 = -A(K_,K + )\ B 21 = -B\ 2 , (j = 1,2). 
5 Quantum Dissipative Systems 

Open systems are systems that interact with other systems and/or their environment. In classical control 
theory, a general framework for the stability of open systems has been developed [551(511 [561 [57]. This classical 
theory abstracts energy concepts and provides fundamental relations for stability in terms of generalized 
energy inequalities. In this section we briefly review dissipation theory for open quantum systems, |23] , In 
Fig. [8] the open quantum system P = (S, L, H) is the plant of interest whose space of variables is denoted 
Ap. The other open quantum system W = (R, w, D) is an external system or the environment, whose space 
of variables is denoted A ex . W is called an exosystem. Moreover, we allow W to vary in a class of such 
exosystems W. 

5.1 Dissipativity, Stability, Passivity, Gain 

In this section we present concepts of dissipativity, stability, passivity and gain. As with the classical case, 
criteria of stability, passivity and gain follow those of dissipativity. 

The following assumption is used in the sequel. 
Assumption Al. The inputs to the composite system P AW are all canonical vacuum fields, cf. Section [3.11 

Let rp{W) be a self-adjoint operator in the composite plant-exosystem space Ap®A ex - rp(W) is usually 
called supply rate. We have the following definition of dissipativity for open quantum systems P. 

Definition 1 (Dissipation, ]23\ Sec. Ill- A]) The plant P is said to be dissipative with supply rate rp{W) 
with respect to a class of exosystems W if there exists a non-negative plant observable V £ Ap, called storage 
function, such that the dissipation inequality 

E [V(t) - V(0)} < [ E [r P (W)(s)} ds (33) 
v v ' Jo t 

stored energy v v 

supplied energy 

holds for all W G W and all t > 0, where Kq is vacuum expectation fJA Chapter 26]. In particular, when 
"=" in ( Eg) holds for all W G W and all t>0,Pis called lossless. 

The combination of Definition [1] and the following property of vacuum expectation |4H Proposition 26.6] 

E T [V(t)} = V(t) + ^ E [g PAW (V(r))]dr, (34) 
yields an infinitesimal version (namely independent of the time variable) of the dissipation inequality ([33|) . 
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Theorem 2 (Dissipation, \23\ Theorem 3.1]) Pertaining to Fig. the plant P is dissipative with a supply 
rate rp(W) with respect to a class of exosy stems W if and only if there exists a non-negative plant observable 
V G Ap such that 

GPAw(V)-r P (W) <0 (35) 

holds for allWeW. 

Given a quantum system P, assume it has indirect and/or direct connections to external systems and/or 
its environment. With slight abuse of notation we still call P an open quantum system. Let W u denote the 
class of all the quantum systems that can be connected to P, directly or indirectly. The following result 
shows that P is lossless with respect to W u - 

Theorem 3 ( \23\ Theorem 3.3]) Pertaining to Fig. for any given storage function Vq, which is a non- 
negative observable in Ap, the open quantum system P is lossless with a supply rate 

r (W) = Gp aW (Vo) (36) 

with respect to W u . 

Theorem [3] shows that any open quantum system is dissipative in some sense. 

Next we study stability of open quantum systems which is characterized in terms of the evolution of 
mean values. 

Definition 2 (Exponential stability, 1231 Sec. III-B]) An open quantum system P is said to be exponentially 
stable if there exists a non-negative observable V G Ap, scalars c > and A > such that 

(V(t)) < e^iV) + - c (37) 

holds for any plant state and all time t > 0. Moreover, if A = 0, then limt_ >00 (y(t)) = 0. 

The combination of Eq. (|34p and Definition [2] gives the following stability result for open quantum 
systems P. 

Theorem 4 (Stability, \23\ Lemma 3.4]) If there exists a nonnegative observable V G Ap, scalars c > 
and A > such that 

g P {v) + cv<\, 

then the open quantum system P is exponentially stable. Moreover, if X = 0, then lim(_ s . 00 (y(t)) = 0. 

In what follows we focus on the series product of P and W with additional direct coupling, Fig. [9j That 
is, the composite system is 

PAW = (P<W)m(-,-,H int ), 

where P = (I,L,H), W = (l,w,0), H int = -i(M^v - v^M) with w,v G A ex and M G A P . With slight 
abuse of notation, we write 

P a W = P < W, 

where W = (I, w, —i(M^v — v^M)). That is, direct coupling is absorbed into the exosystem W. 
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Figure 9: Series product plus direct coupling 
Let V G Ap be a non-negative observable. Assume that S = I. By Eq. ([14 



GpawW = Q p {V) +C w (V) + [wi v^]Z + Z^ 



w 

V 



+ [V,v^]M - M^[V,v] 



(38) 



where Z 



V. 



L 

M 



For fixed M E _4p, define a class of exosystems 

u — v^M) : w, v commute with Ap}. 
Then we have the following definition of passivity for the system P in Fig. [9l 



(39) 



Definition 3 (Passivity, \23\ Sec. III-C]) Given M £ Ap, the plant P = (I,L,H) is said to be passive 
with respect to the class of exosystems Wi in $39\) if it is dissipative with the supply rate 



r P (W) = -N^N + [w ] v ] ]Z + Z ] 



w 

V 



+ A 



(40) 



for some non-negative real number X, and N,Z G Ap. P is said to be strictly passive if N^N is strictly 
positive. 

The combination of Definition [3l Theorem [2] and Eq. ()38[) gives the following passivity result. 

Theorem 5 (Positive Real Lemma, I2R Theorem 3.6]) A plant P = (I, L, H) is passive with respect to the 
class of exosystems Wi in I^Pj) if and only if there exists a non-negative observable V £ Ap, an operator 
N G Ap, and a non-negative real number A such that 



Qp(V) + N^N - A < 

Z -- 

As with the classical case, strict passivity implies stability. 





L 




v, 




M 





(41) 
(42) 



Theorem 6 The open quantum system P is exponentially stable if it is strictly passive with respect to the 
exosystem W = (1, 0, 0). 

The bounded real lemma is used to determine L 2 gain of an open system, and in conjunction with the 
small gain theorem, can be used for robust stability analysis and design. In what follows we discuss L 2 gain 
of open quantum systems. 

Define a class of exosystems 



W2 = {W = (I, w, 0) : w commutes with Ap}. 
(Note that in this case there is no direct coupling.) 



(43) 
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Definition 4 (1? gain, \2S\ Sec. III-CJ) The plant P = (I,L,H) is said to have I? gain g > with respect 
to the class of exosystems W2 in ft43\ ) if it is dissipative with the supply rate 

r P {W) = g 2 w ] w-(N + Zw)\N + Zw) + \ (44) 

for some non-negative real number \, and N, Z £ Ap. 

The combination of Definition [3] and Theorem [2] gives the following result. 

Theorem 7 (Bounded Real Lemma, \2S\ Theorem 3.7]) A plant P = (I,L,H) has L 2 gain g > with 
respect to W2 if and only if there exists a non-negative plant variable V £ Ap, an operator N € Ap, and a 
non-negative real number A such that 

r = g 2 - Z^Z > (45) 

and 

Gp{V) + N^N - w^Tw + w\[V, L] + Z*N) + ([V, L] + Z^Nfw - A < (46) 
for all w G A ex . IfT~ l exists, then plant P = (I, L, H) has gain g > with respect to W2 if and only if 

Gp(V) + N^N + {[V, L] + ^J\O t r~ 1 ([V J L] + Z ] N) - A < (47) 

for all w 6 A ex - In the latter case, the plant P is strictly bounded real. 

5.2 Example 

The following example illustrates the above results for stability, passivity, and L 2 gain. Consider a two-level 
system P and an exosystem W of the form 

P = (l )X /7<7 + ,|a 2 ), W = (l,w,0), 

where w commutes with Ap. Assume there is no direct coupling between P and W. Choose a storage 
function Vq = \{I — a z ) = a-(T+ and a supply rate 

r P (W) = g PAW (V ) = Qp{Vo) + w*Z + Z*w, (48) 

where Z = [Vb, ^/jo- + ] = —^/ja + . Clearly, Qp(Vo) = — 7V0. As a result, Eq. (138]) becomes 

rp(W) = — 7V0 — y/j(w*a + + a-w) 

= —(y/j<T+ + w)*(^/jCT + + w) + W*W. 

Choose N = y/jcr + and Z = —N, by Theorem [5] we see that the system is passive. Choose = ■ v /7<t+ 
and Z = 1, by Theorem [7] we find that the system has L 2 gain 1. Finally, when W = (1,0,0), rp(W) = 
rp(I) = — 7V0, then by Theorem [3J the system P is exponentially stable. 

6 Linear Quantum Systems 

Linear quantum systems are those for which certain conjugate operators evolve linearly, the optical cavity 
being a basic example, cf. Section 13.3.11 Linear systems have the advantage that they are much more 
computationally tractable than general nonlinear systems, and indeed, powerful methods from linear algebra 
may be exploited. 
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Figure 10: General model 

6.1 General Model 

Open linear quantum systems discussed in this paper are open quantum harmonic oscillators with direct 
and indirect couplings to other quantum systems and/or external fields. In this section we present a general 
model for an open linear quantum system G, Fig. [TU1 based on the ingredients discussed in the previous 
sections. Here, G is an open quantum system with parametrization (I, L, Hq), where L = C-a + C+a& with 
C_,C+. being constant complex- valued matrices. The internal Hamiltonian Hq is that given in Eq. (0). 
Moreover, G is allowed to coupled directly to another (independent) quantum system W d via an interaction 
Hamiltonian 

H mt = - (a^v + ^Ha) , (49) 

where E = A(iK-, iK+). Our interest is in the influence of external systems/fields on the given system G. 
The performance characteristics of interest are encoded in a performance variabl^ z. 

Building upon the discussions in previous sections, the equations for G (including direct coupling, indirect 
coupling and performance variable) are 

S(t) = Aa(t) + B d v{t) + B f w(t) + B f b(t), a(0) = a, (50) 
bout(t) = C f a(t) + w(t) + b(t), (51) 
z{t) = C p a(t) + D pd v(t) + D pf w(t), (52) 

The complex matrices in (150p and (I5ip are given by 

A = -ic}C / -A(iO_,iO + ), B d = -A(K_,K+)\ (53) 

C f = A{C- ) C+), Bf = -C). (54) 

The matrices A and Bf are specified by the parameters Q±, and C±. In equations (|50p and (|5ip . b(t) and 
b ut(t) are respectively, the input and output fields for G. The term v in (|50l) is an exogenous quantity 
associated with W d with which G is directly coupled via the interaction Hamiltonian Hint > cf . Section 14.51 
The term w in (|50p is another exogenous quantity associated with another (independent) system Wf 
with which G is indirectly coupled through a series product. Wf may be a quantum system of the form 
(l,w, 0) where w is an operator on some Fock space, it can also denote modulation so that w coherent drive 
modulates the vacuum field b cf. E in Example 2 of Section 14.51 Because of the assumed independence, 
w and v commute with the mode operators dj,a* for G. While v and w are arbitrary external variables, 

2 A performance variable is chosen to capture control performance, such as an error quantity, and so may involve external 
variables, like a reference signal. Performance variables need not have anything to do with the output quantities associated 
with direct or indirect couplings to other systems, cf. Sections [7] and [8] 
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the time evolutions v(t) and w(t) (when it is an operator) are determined by the evolution of the overall 
composite system. The matrices C p , D p d and D p f specify the performance variable z. In brief, system 
G is specified by the parameters G = (£l±,C±, K±,C P , D p d, D p f). Of these, Q±, C± and K± are physical 
parameters. 

In particular, when all the plus terms are zero, namely C+ = 0, n + = 0, K + = 0, all matrices 
A, Bf, Bd,Cf are block diagonal, system ([50]) - ([5T]) is equivalent to 



a(t) = 
bout(t) = C-a(t) + b(t) 



(in. + ^ciCL)a(f) - K ] _v{t) - C ] _w(t) - Ch(t), a(0) = a, 



(55) 
(56) 



cf. optical cavity (|19 p - (|20p . It can be readily shown that system (|55D - ([56l) is passive. Passive systems have 
been studied in, eg., [281 E3 EH ESI M\ . 



6.2 Physical Realizability 

It can be readily verified that the following relations for system matrices (|53p -([54|) hold 

Jn,A + J n + C^J m Cf = 0, 



B 



f 



-G b 



B d = -A(K_,K + f. 
Equation (|57p characterizes preservation of the canonical commutation relations, namely 

[aj(t),a* k {t)} = [aj,a* k ] = (J n ) jk , Vi > 0, (j, k = 1, . . . , n). 



(57) 
(58) 
(59) 

(60) 



Eq. (I58p reflects the input and output relation, while Eq. (I59|) is for direct coupling. 

The relations (|57p -(|59p are called physical realizability relations, which generalize results in \12\ Theorem 
3.4], [56], [291 Theorem 5.1], Theorem 3]. These conditions guarantee that the equations correspond to 
a physical system. 



6.3 Quadrature Representation 

So far, annihilation- creation representation has been used to represent linear quantum systems in terms of 
the notation a = [a T a^] T , the resulting matrices are complex-valued matrices. In this section we introduce 
an alternative representation, the so-called quadrature representation, which leads to equations with real- 
valued matrices. 



Define the unitary matrix 



and the vector of self-adjoint operators 



A 



1 

71 



i i 

-il il 



Q 
P 



by the relation 



5 = Aa. 



(61) 



(62) 



(63) 



The vector q = -^[1 I]a is known as the real quadrature, while p = "^[ — H *^]& is called the imaginary or 
phase quadrature [53] . 
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Similarly define unitary matrices Af, A d and A p of suitable dimension, of the form (|6ip . and define 
quadrature vectors 

b = Afb, b out = Afb out , w = Afw, v = A d v, z = A p z. 
Then in quadrature form G is in the form 



a(t) = Aa(t) + B d v(t) + B f w(t) + B f b(t), 5(0) = a, 

(t) = C f a(t) + w(t) + b(t), 
z(t) = C p a(t) + D pd v(t)+D pf w(t), 



(64) 



where A = AAA" 1 ", B d = AB d A d , B f = AB f A\, C f = A f C f A\ C p = A p C p A\ D pd = A p D pd A\, D pf = 
A p D p fA^. Note that all entries of the matrices in this representation are real. 

6.4 Series Products for Linear Quantum Systems 

Assume both G\ and G2 in Fig. H]are linear, in this section we present the explicit from of G = G2<G±. 

For ease of presentation we assume both Gi and G2 are passive with parametrization Gj = (I, C_ 0), 
(j = 1,2). Therefore 



C®aW(t) + b®(t), (j = 1,2). 



According to Eq. (j29|) . the composite linear quantum system G is 



G 



I, [C W C m ] 



a< 2 > 



2i 



[(«W)t (a«)t] 





(C (2) )tc (l) 



, (C (D )tC (2) 



By Eqs. (|55p-(|56p. G is in the form of 



dW(t) " 




" |( C (i))tc^) 




" aW(i) " 




(C (1) )t 


a( 2 )(t) 




(c (2) )tc7 (i) (G P))tc7i 2 ) 




a( 2 )(t) 




-(ci 2) )t 



a( 2 ) 



6W(t) 



(65) 
(66) 

(67) 



[c a) c (2)] 



a«(t) 
a( 2 )(i) 



+ &«(*). 



(68) 
(69) 



If we identify ft( 2 )(t) with 6^(t) in Eqs. (|55]l -([5B jl . then Eqs. (|55]l -((55 ]l give rise to system ((SHI)- HMD to °- 



This fact is useful in forming closed-loop coherent feedback control systems, cf., Section 1871 



7 Performance Specifications for Linear Quantum Systems 

In Section [5] we have established criteria for stability, passivity, and 1? gain for general quantum dissipative 
systems studied in Section 13.21 These criteria are expressed in terms of operators. In this section we 
specialize those results to linear quantum systems introduced in Section (SI It can be seen that for linear 
quantum systems such criteria can be expressed in terms of constant matrices. 
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7.1 Stability, Passivity, Gain 



Perhaps the most basic performance characteristic is stability. For system G of open quantum harmonic 
oscillators presented in Section 16.11 stability may be evaluated in terms of the behavior of the number of 
quanta (e.g. photons) stored in the system, N = a} a = ]Cj=i a *j a j- We introduce the following definition of 
stability. 



Definition 5 (Stability, \6(^ Sec. Ill- A]) Let w = and v = in 150\) . that is there is no energy input 
to system G. We say that G is (i) exponentially stable if there exist scalars cq > 0,ci > 0, and c 2 > 
such that (NO)) < coe~ Cl '(N) + c 2 ; (ii) marginally stable if there exist scalars c\ > and c 2 > such that 
(N(i)) < ci(N) + c 2 t; and (Hi) exponentially unstable if there exists an initial system state and real numbers 
c > 0, ci > and c 2 such that (N(i)) > c e Cl *(N) + c 2 . 



For example, for the closed optical cavity in Section l3.3.H Fig. [U a(t) = exp(—iu)t)a, and a*(t)a(t) = a* a 
for all t, which means that G is marginally stable but not exponentially stable — it oscillates — hence the name 
"oscillator". However, an open cavity (1, ^/na, cja*a) (Fig. [2]) is exponentially stable, a damped oscillator. 

The number operator N = a) a, whose mean value is the total number of quanta, is a natural Lyapunov 
function for G, and is directly related to the energy of the system. However we find it more convenient 
to use storage functions of the form V = ^a^Pa for non-negative Hermitian matrices P. For such storage 
functions, the generator function (|14p becomes 



Qa{V) = h\A^P + PA)a. 



(70) 



With this simple yet important observation, the results in Section can be specialized to linear quantum 
systems. 

Define the matrix F by 



Fdt = {dB*{t)dB T {t)) ri 



m 

I m 



dt. 



(71) 



The following result is a simple criterion for stability of linear quantum system G, which is a linear 
version of Theorems H] and [6l 

Theorem 8 (Stability, JUB, Theorem 1]) If there exist constant matrices P > and Q > cP for a scalar 
c > such that 

A ] P + PA + Q<0, (72) 

then inequality 



a!(t)Pa(t)) < e 



a l Pa) + 



Yc 



(73) 



holds, where A = tv[B\PB f F] with F given by j fTip . If also P > al (a > 0), then (a t (t)a(t)} < 



i e -ct 



a) Pa) + t£-. In this case, G is exponentially stable. 



In a similar way, by choosing linear versions of supply rate functions positive real lemma and bounded 
real lemma can be established for linear quantum systems. 

In order to simplify the notation we write u = [w T v T ] T for the doubled-up vector of external variables, 
and define accordingly 

" I 



B := [B f B d ] 



0/0 
0/00 
/ 



(74) 
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where dimensions of identity matrices are implicitly assumed to be conformal to those of v and w. 
Define a supply rate 

r(a, u) = -(— a'Qa + ii) z + z*u). 
The we have the following positive real lemma. 



(75) 



Theorem 9 (Positive Real Lemma, \66[ Theorem 3]) The system G with performance variable z = C p a is 
passive if and only if there exist non-negative definite Hermitian matrices P and Q such that 



PA + A^P + Q PB-Cl 
B^P - C„ 



< 0. 



(76) 



Moreover, A = tr[B j f PB f F]. 
Remark 3 When 

P = H - 



n_ n + 
n* n* 



, V = -tfPa, L = C-a + C+a*, M = K^a + K + a#, 



the operator Z in Eq. (|42|) satisfies Z = z = C p a. That is, Theorem is a special case of Theorem 

In what follows we discuss 1? gain of linear quantum systems. Denote D p = [D p f D p( j\, the performance 
variable can be rewritten as z = C p a + D p u. Define a supply rate 



r(a, u) = - - [pz - g 2 u^u), 
where g > is a real gain parameter. The we have the following results. 



(77) 



< 0. 



(78) 



Theorem 10 (Bounded Real Lemma, \bb\ Theorem 4j) The system G with performance variable z = C p a + 
D p u is bounded real with finite L 2 gain less than g if and only if there exists a non-negative Hermitian 
matrix P such that 

PA + A^P + clC p PB + C ] P D P 
B^P + D ] P C P D ] p D p - g 2 L 

Moreover, A = tv[B^PB f F]. 

Theorem 11 (Strict Bounded Real Lemma, \6b\ Theorem 5]) The following statements are equivalent. 
i) The quantum system G defined in i5(J\) -[o ^) is strictly bounded real with disturbance attenuation g. 



ii) A is stable and C p (sL - A) 1 B + D t 



< 9- 



iii) g 2 L — D P D P > and there exists a Hermitian matrix P\ > satisfying inequality 



A^P 1 +P 1 A P X B C\ 



B^P X 



-gi 4 
-gi 



< o. 



(79) 
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iv) g 2 I — DpD p > and there exists a Hermitian matrix P 2 > satisfying the algebraic Riccati equation 

A^P 2 + P 2 A+( P 2 B + (AD. 



2 -t- r 2 A f \tid f 

x - £>{!>,) ~* (fitpt + Dtc p ) 
= 

with A + BB^P2 being Hurwitz. 
Furthermore, if these statements hold, then P\ < P2 . 

7.2 LQG Performance 

In this section a quantum LQG cost function is first denned in the annihilation-creation form, and whose 
evaluation is connected to a Lyapunov equation in the complex domain. After that the real domain case is 
presented. More discussions can be found in, e.g., [371 [66], 
Consider the following stable linear quantum system 

da{t) = Ad{t)dt + B f dB(t) (80) 

where B(t) is a quantum Wiener process introduced in Section [3.11 Given a performance variable z{t) = 
C p d(t), along the line of [37], the infinite-horizon LQG cost is 

3oc := Um -J- 1' \(z\t)z{t)+z T {t)z#{t))dt (81) 

t/-K» tf J 2 \ I 

= Tr {C P P LQG CI} , 

where the constant Hermitian matrix Plqg ^ satisfies the following Lyapunov equation 

APlqg + PlqgA^ + \b s b) = 0. (82) 

In quadrature form, given a stable linear quantum system 

da(t) = Aa(t)dt + B f dB(t) (83) 

with performance variable z(t) = C p a(t). Assume that the constant real matrix Plqg > is the (unique) 
solution to the following Lyapunov equation in the real domain 

APlqg + Plqg A T + B f Bj = 0. (84) 

Then 

Zoo = ^[c p P L qgCI). (85) 
8 Coherent Feedback Control 

We have discussed interconnections of quantum systems (Section U]), open linear quantum systems (Section 
[6|), and their performance specifications (Section [7J). We are now in a position to study synthesis of open 
linear quantum systems; that is, how to connect a plant of interest to another system (namely controller) 
so as to achieve pre-specified control performance. 
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Figure 11: Coherent feedback control arrangement 



8.1 Closed-Loop Plant-Controller System 

In Figure [TTl P is the plant to be controlled, K is the controller to be designed. Clearly, this feedback 
system involves both direct and indirect couplings between P and K. 

The plant P is described by a system of quantum stochastic differential equations (QSDEs) 



S(t) = Aa{t)+B 12 a K (t) + B v b v {t) + B f w{t) 

+B f b(t) + B u u(t), 8(0) = a, 
y(t) = Ca{t) + D v b v {t)+D f w(t) + Dfb(t). 



(86) 



The inputs w(t) and b(t) are defined in Section [(J y(t) is a selection of output field channels from the plant. 
b v (t) is a vector of additional quantum white noises; u(t) is a quantum field signal from the to-be-designed 
controller K, hence it is a vector of physical variables. The term Bi 2 aK(t) is due to direct coupling between 
P and K. 

The fully quantum controller K is a linear quantum system of the forrrj§ 



CLK(t) 



A K a K {t) + B 2l a{t) + B K y(t) + B KX b VK1 (t) 
+B K2 b VK2 (t), a K (0)=a K , 
u{t) = C K a K {t) + b VK1 (t). 



(87) 



This structure allows for direct coupling and indirect coupling between the plant P and the controller K. 
Here, b VK1 {t) and b VK2 (t) are independent quantum white noises, and u(t) is the field output of the controller 
corresponding to b VK1 {t). Finally the terms B^axit) and B 2 ia(t) are due to the direct coupling between 
the plant and controller in terms of an interaction Hamiltonian 



Hm = ^ {a ] r)a K + 8^H8 



where S = A(iK-, iK+) for complex matrices K_ and K + of suitable dimensions, cf. Section [4.31 

The controller matrices K + , (or B± 2 , B 2 \) for direct coupling, and Ak, Bk, Ck, Bki, Bk2 for indirect 
coupling are to be found to optimize performance criteria defined in terms of the closed-loop performance 
variable 

a(t) 



m = [c P d u c k ] 



a K (t) 



+ D pf w(t). 



(89) 



We assume that all the variables and matrices of the plant and the controller have compatible dimension, but we don't 
bother to specify them explicitly. 
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Because standard matrix algorithms will be used in H°° synthesis and LQG synthesis in later sections, 
we resort to quadrature representation discussed in Section 



the quadrature counterparts of a, &k, w, b, b v , z, {3 U , y, b 



V K1 > U V K2 



Let a, ax, w, b, b v , u, z, y, b VK1 , b VK2 be 
respectively. Define 



A B u Cr 
BrC A k 



+ 3, B, 



d 



B K Df 



G c i 



Bf B v B u 
BrDj BkD v Bki Bk2 



C c l 



Cp D u Ck 



D c i — D p f, 



where S = [0 B\2] B21 0] satisfies B21 = QBj^Q. Then the closed-loop system in the quadrature represen- 
tation is given by 



a(t) 

a K (t) 



5(f) 
a Kit) 

ait) 
a K it) 



+ B d w(t) + G c 



+ D d w(t). 



bit) 

Kit) 

b VK1 it) 
bv K2 it) 



(90) 



(91) 



8.2 H°° Control 

As in the classical case, the bounded real lemmas stated in Section [7TT1 can be used for H°° controller synthesis 
of open linear quantum systems. It is shown in [22] that for open linear quantum systems H°° control 
performance and physical realizability condition of controllers can be treated separately. Adding direct 
coupling between plants and controllers complicates H°° controller synthesis. Nonetheless, the separation of 
H°° control performance and physical realizability condition still holds. This is a unique feature of quantum 
H°° controller synthesis: To guarantee physical realizability, vacuum noise is added, while such noise does 
not affect H°° control performance [22]. 



8.2.1 LMI Formulation 

In this section we present a general formulation using LMIs for H°° synthesis of open linear quantum 
stochastic systems. 

According to the strict bounded real lemma (Theorem II ip . the closed- loop system (|90p -(|9ip is internally 
stable and strictly bounded real (from w to z) with disturbance attenuation g if and only if there is a real 
symmetric matrix V such that 



V > 



AJ,V + VA cl VB cl CT, 



-gi 



-gi 



< 0. 



(92) 
(93) 



The H°° controller synthesis is to find indirect coupling parameters Ak, Bk, Ck and direct coupling 
parameters S such that Eqs. (f92|) - ([93|) hold. 
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V 



Partition V and its inverse V 1 to be 

Y N 

N T * 

Define matrices 

X I 

M T 



V 



ni 



n. 



X M 

M T * 

I Y 

N T 



And also define a change of variables 

A = N{A K M T + B K CX) + Y(B U C K M T + AX), 
B = NB K , 



C 

n 



CkM 1 , 



With these notations, ([92]) and ([93]) hold if and only if the following inequalities hold. 

< 0, 



X I 
I Y 



* 


* 


* 


* 


-gi 


* 


Del 


-gi 



+ 



AX + Xi T + B U C+(B U C) T 
A+A T 

CpX+Z? u C 

A + A T 
A T Y + YA + BC + (BC) T 
{YB f + Kb f ) T 

' B 12 M T + {B 12 M T ) T {NB 21 X) T + {YB 12 M T f 
NB 21 X + YB 12 M T NB 21 + (A^S 21 ) T 











< 0. 



(94) 



(95) 



(96) 



If (|95p and (|96p are simultaneously soluble, according to Eq. (|94p . the following matrices can be obtained. 

B K = iV^B, (97) 

C K = c(M T y\ 

A K = N- 1 (k-NB K CX-Y(B u C K M T + AX))M- T , 
E = V~ l (a\ T ) flUi 1 . (98) 

Unfortunately, notice that there are such terms as NB 2 \~K and YB\ 2 M T in inequality (|96p . which induce 
nonlinearity. The above analysis shows it is hard to directly utilize LMI techniques to do controller design 
when direct coupling is involved. 
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8.2.2 Multi-step Optimization 

In this section, we attempt to circumvent the above difficulty by proposing a multi-step optimization pro- 
cedure which is formulated as follows: 

Initialization. Set B\2 = and B21 = 0. 

Step 1. Solve linear matrix inequalities (j95|) and ([96]) for parameters A, B, C, X, Y and disturbance 
gain g, then choose matrices M and ./V satisfying MN^ 1 = I — XY. 

Step 2. Pertaining to Step 1. Solve inequality (1961) for direct coupling parameters B\2,B 2 \ and 
disturbance gain g. 

Step 3. Fix B12 and B 2 \ obtained in Step 2 and M and iV in Step 1, go to Step 1. 

After the above iterative procedure is complete, use the values B\2,B2i, Ak, Bk,Ck obtained to find 
Bri,Bk2 to ensure physical realizability of the controller. A complete procedure of finding matrices 
Bki,B K2 is given in [221 Sec. V-D]. 

Remark 4 Steps 1 and 2 are standard LMI problems which can be solved efficiently using the Matlab 
LMI toolbox. However, there is some delicate issue in Step 3. Assume that B\2 and B21 have been obtained 
in Step 2. According to the second item in (|96j) . constant matrices M and N must be specified in order to 
render (|96p linear in parameters A, B, C, X, Y, and disturbance gain g. In Step 3, M and N obtained 
in Step 1 is used. Unfortunately, this choice of M and N sometimes may generate a controller whose 
parameters are ill-conditioned. Due to this reason, M and N in Step 3 might have to be chosen carefully 
to produce a physically meaningful controller. This fact is illuminated by an example in [661 Sec. IV-C6]. 

Finally we discuss robustness briefly. It is demonstrated in [66] that direct coupling may improve 
robustness of closed-loop quantum feedback systems. For instance, for the example studied in [22J Sec. 
VILA], using coupling coefficients n\ = 2.6, K2 = K3 = 0.2, we implement Step 1 of the above multi-step 
optimization procedure to design an indirect coupling, and obtain closed-loop L 2 gain 0.0487. We implement 
Step 2 to design direct coupling and obtain an L 2 gain of 0.0498. This is a bit worse than the previous one, 
however the difference is quite small. Now we assume there is uncertainty in the coupling coefficient Ki, 
say the actual value of K\ is 1.3. In this case, the L 2 gain of the closed-loop with indirect coupling becomes 
0.1702, which is a significant performance degradation. However, the L 2 gain of the closed-loop with both 
direct and indirect couplings is 0.0595, which is still close to the original 0.0498. 

8.3 LQG Control 

In this section we study the problem of coherent quantum LQG control by means of both direct and indirect 
couplings. In contrast to the coherent quantum H°° controller synthesis presented in Section 18.21 the nice 
property of separation of control and physical realizability does not hold any more. This is evident as 
LQG control concerns the influence of quantum white noise on the plant, the addition of quantum noise 
that guarantees the physical realizability of the to-be-designed controller affects the overall LQG control 
performance. 

In the following we just give a brief formulation of the coherent quantum LQG control problem. In-depth 
discussions can be found in [371 ESJ E3j [M] • 

We make the following assumption. 
Assumption A2. There are no quantum signal w{t) and noise input b v (t) in the quantum plant P in (|86p . 

Following the development in Section 17.21 the LQG control objective is to design a controller (|8T|) 
such that the performance index 3oo = Tr |c' c /PxqgC'^| is minimized, subject to equation (|84p and the 
quadrature counterpart of the physical realizability condition (|57l) - (|59l) . 
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As yet, quantum LQG coherent feedback is still an outstanding problem, there are no analytic solutions. 
In |37] an indirect coupling is designed to address the coherent quantum LQG control problem, where a 
numerical procedure based on semidefinite programming is proposed to design the indirect coupling. In 
order to design both direct and indirect couplings. In [66\ Sec. IV-D] a multi-step optimization algorithm 
is developed to incorporate direct coupling into numerical design procedures. Some theoretical insights into 
the structure of quantum LQG coherent feedback control have been provided in [63^ I64j . 

9 Network Synthesis 

A linear quantum controller, obtained from either coherent H°° or LQG control synthesis, is in the form of 
a set of linear quantum stochastic differential equations. Network synthesis theory is concerned with how to 
physically implement such controllers by means of physical devices like optical instruments. This problem 
has been addressed in [38\ \39\ [4"0"1 [4"2] . The general result is: A general linear quantum dynamical system 
can be (approximately) physically implemented by linear and nonlinear quantum optical elements such as 
optical cavities, parametric oscillators, beam splitters, and phase shifters. 

Lately, the Mabuchi group at Standard [50] has developed a Quantum Hardware Description Language 
(QHDL) to facilitate the analysis and synthesis of quantum feedback networks described in this survey. 
As a subset of the standard Very High Speed Integrated Circuit (VHSIC) Hardware Description Language 
(VHDL), QHDL provides high-level modular representations of quantum feedback networks. This user- 
friendly interface will be helpful to the fabrication of complex photonic circuits. 

10 Conclusions 

In this survey we have presented a brief look at recent results concerning quantum feedback networks and 
control. On the basis of this model interconnection structures of quantum systems have been presented. 
Fundamental characteristics of quantum systems such as stability, passivity, and 1? gain have been described. 
It turns out that for linear quantum systems these fundamental characteristics have very explicit forms. The 
problem of coherent H°° control and coherent LQG control have been discussed. 
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